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Last years many interesting effects related to incoherent MAR have been experimentally found, but only 
few of them were theoretically explained. It was shown, for example, that if the voltage at the edges of a 
linear array is V then subgarmonic structures in the current -voltage characteristics appear not only at usual 
for nonstationary Josephson effect positions, V„ — 2A/n, where n is integer, but also at voltages other than 
V„. A step towards description of electron transport in a dirty array of SNS junctions is done in this letter. 
It is shown that subgarmonic structures may indeed appear at "unusual" voltages. 
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Important role plays Andreev reflection mechanism 
in the subgap charge transfer through a normal metal 
(N) - superconductor (S) junction p]. When an electron 
quasiparticle in a normal metal with the energy below 
the superconduting gap reflects from the interface of the 
superconductor into a hole. Cooper pair transfers into 
the superconductor. If the normal metal is surrounded 
by superconductors, so we have a SNS junction, a num- 
ber of Andreev reflections appear at the NS interfaces. 
In equilibrium this leads to Andreev quasiparticle levels 
in the normal metal that carry considerable part of the 
Josephson current; out of the equilibrium, when super- 
conductors are voltage biased, quasiparticles Andreev 
reflect about 2A/eV times transferring large quanta of 
charge {r-^ 2e • [2A/eV]) from one superconductor to the 
other. This effect is called Multiple Andreev Reflec- 
tion (MAR). If the voltage is near Vn = 2A/n, where 
n = 1, 2, . . ., so-called subgap features in current voltage 
characteristics appear. Then large contribution to the 
current give quasiparticles that go from the gap edge 
of one superconductor to the gap edge of the other su- 
perconductor (after MAR in the normal region); bulk 
superconductor DoS is large at the gap and this is the 
reason of subgap features. 

This letter is devoted to investigation of electron 
transport in arrays of dirty superconductiong meso- 
scopic SNS (SFS) junctions. I assume that normal 
parts of the junctions are "long" . It implies that min- 
imum distance do between adjacent superconductors is 
much larger than the characteristic scale of anomalous 
green function (Cooper pair wave function) decay 
from a superconductor in the normal metal. If the 
diffusion coefficient of the normal metal is Dn than 



^ \/Dn/T ^ do, where T is the temperature. En- 
ergy relaxation is not included in the calculations. So it 
is also assumed that the array is not large, its charac- 
teristic length does not exceed the characteristic length 
scale of quasiparticle energy relaxation in the normal 
metal and in superconductors. [Quasiparticle energy 
relaxation was "taken into account" in Ref. [1], where 
incoherent MAR in long SNS junction was discussed, 
by small imaginary part supplied to the energy in re- 
tarded and advanced greens functions; but coUisional 
integrals were not taken into account in kinetic equa- 
tions. It is not clear why this procedure is correct. I do 
not follow this way here.] Conditions listed above also 
mean that in equillibrium the proximity (Josephson) ef- 
fect between the superconducors is suppressed. Joseph- 
son current, for example, is exponentially small with 
do/^N- It is known that out of the equilibrium when 
there is a finite bias between the superconducors the 
proximity effect restores in some sense: subgap features 
appear due to MAR in current-voltage characteristics. 
MAR in long Josephson junctions usually is referred to 
as 'incoherent" since there is no contribution to electron 
transport from effects related to interference of quasi- 
particle wave functions in normal metals (contrary to 
Josephson effect in "short" superconducting junctions). 
Last years many interesting effects related to incoherent 
MAR have been experimentally found, but only few of 
them were theoretically explained. It was shown, for ex- 
ample, that if the voltage at the edges of a linear array 
is V then subgarmonic structures in the current -voltage 
characteristics appear not only at usual for nonstation- 
ary Josephson effect positions, Vn = 2A/n, where n is 
integer, but also at voltages other than Vn, see Ref. [5] 
and refs. therein. A step towards description of elec- 
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Fig.l. a) An array of SNS junctions, like in experi- 
ments [3 . b) A sketch of the normal layer connection to 
the superconductors in experimental SNS arrays. Black 
area is insulating. 



tron transport in a dirty array of SNS junctions is done 
in this letter. It is shown that subgarmonic structures 
may indeed appear at "unusual" voltages. Investigation 
of electron transport is based on Usadel equations: 

[H,G]=iDV-J, J-GVG, 6"^ = !, (1) 
ff- 1(1(7,9* -60 + A), (2) 

I -A* 0/ 
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where 



G = 



R K' 
, A, 



(3) 
(4) 

(5) 



Here i?, K and A denote retarded, Keldysh and ad- 
vanced quasiclassical green functions. The hat reminds 
that greens function are in turn matrices in Nambu 
space, a are Pauli matrices that act in Nambu space; 
I? is a diffusion coefficient that equals Djq in a normal 
metal and Ds in a superconductor; A is the supercon- 
ducting gap; (j) is electrical potential; crjv is the conduc- 
tivity of a normal metal. The boundary conditions for 
the Usadel equations at NS interfaces are: 



cr«2 J2 ■ n2 = ^^[62, Gi]. 



(6) 



where 1,2 label interface sides, 6*12 is the surface con- 
ductance and n2 is the unit normal to the interface 
pointing to the second half-space. 

The problem is to calculate the current in a long SNS 
array using Eqs. ((!])- ([5]). Solving these equations directly 




Fig.2. The simplest array of SNS junctions: SNSNS. 



is a hard task because they are nonlinear, nonuniform 
and essentially time dependent (relative phases of su- 
perconductors rotate with biases). So the main task is 
developing an approach for the problem in hand that 
allows to perform significant part of transport calcula- 
tions analytically and that is applicable in rather wide 
range of system parameters. There is no universal ap- 
proach that helps to solve Usadel equations analytically, 
so any analytical method of Usadel equations solution 
is usually specific to the given class of the physical sys- 
tems. 

Normal layers in experimental SNS arrays [2^ con- 
nect with superconductors like it is shown in FiglTJ SNS 
junctions of this type are usually referred to as "weak- 
links" [3 . Boundary conditions, Eq.®, can be simpli- 
fied in this case: retarded and advanced Greens func- 
tions at superconducting sides of NS boundaries can be 
substituted by Greens function from the bulk of the su- 
perconductors. These "rigid" boundary conditions ap- 
proximation are reasonable because a) the magnitude 
of the current is much smaller than the critical cur- 
rent of the superconductor (this is assumed), b) the 
current entering the superconductor from narrow nor- 
mal metal wire [the width f = yjDs/Tc, where Tc is 
the critical temperature of the superconductor] spreads 
nearly at the NS interface over the whole superconduc- 
tor. There are also other cases when rigid boundary 
conditions are correct, for example, if the NS bound- 
ary has small transparency due to an insulator layer or 
other reasons. 

The Keldysh greens function has the following gen- 
eral parametrization: K — R o f ^ f o A, where / 
is the distribution function. It was shown in Ref. [1] 
that the current in a long SNS junction can be found 
from investigation of the current distribution in an effec- 
tive network where the role of voltages play distribution 
functions made from components of /, the role of resis- 
tances play NS resistances renormalized by proximity 
effect and normal layer resistances. This idea can be 
applied to an SNS array. 

It is convenient to write / — lf+ + azf-- Phases 
of the superconducting order parameters rotate with 
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the biases so one could write in general /(e, e') — 
E«={„,;i=i.2...} /"(e - - e' + niV,), where n, is in- 
teger, Vi is the bias at the z's superconductor. (Similar 
consideration apply for retarded and advanced greens 
functions.) However due to the absence of Josephson- 
type interference effects nonzero are only that /„ which 
provide time independent (dissipative) component of 
the current as it is well explained in Ref. i4j... One can 
write boundary conditions Eq. ^ at an NS interface for 

fn 9-8 



(7) 
(8) 



Here the label 2 corresponds to the normal metal. For 
example, 



D+ = ^Tr(l-RA], 
D_ = i Tr ( 1 - Ra.Aa, 



(9) 
(10) 



Definitions of G± are similar and can be found, e.g., in 
Ref. [4 . 

It is useful to go from f± to ngj-zj) 41 that are related 
as follows; J^. = 1 — (rig + rih), /- = (fih — Ue). Then 
the boundary conditions, Eqs.Q, can be written as 



/g = GT(nW-n(2))~ 

-G^((n(2)-4^))-(nW"4'^)), 
/. = GT(n«-4^^)+ 

+G^((n(2)-ni^))-(nW-nW)). 



(11) 
(12) 



Here Gt ^ G+, Ga = (G_ - G+)/2 and = 
—o'N ^iD+d{ne + Ufi) ± D^d{ne — nh))- If the supercon- 
ducting bank labelled by index "1" is in equilibrium, so 
rie^ — n^j^^ — np, then we arrive at boundary conditions 
written in Eq.(21) of Ref. gj. If one writes Eas. (fTT]) -(ITa 
for NS and SN interfaces of a superconducting island 
and uses conditions of electron and heat currents con- 
servation in the island then quasiparticle distribution 
functions corresponding to superconducting islands can 
be excluded from boundary conditions: 



le 
h 



A, 
A, 



,(2) 



(n(^)-nW) + A_(4^'-nW). 



(13) 
(14) 



It is worth noting that Eqs. (IT5)) - p^ are derived for a 
linear array of SNS junctions; generalization of these 
equations is straightforward for more complicated ar- 
rays. In the last pair of equations labels 1 and 2 corre- 




Fig.3. a) An illustration of the boundary conditions 
Eas. H13|) - ()12|) in terms of a circuit is given in Fig[3] 
Electron and hole currents entering the left side of the 
pyramid flow in one normal layer, the right currents 
flow in the other normal layer. —A- resistance de- 
scribes Andreev reflection, A+ — quasiparticle normal 
transmission through the superconductor and A_ is An- 
dreev transmission characteristics, b) An illustration 
of the boundary conditions Eas. (fT5)) - (fT^ (developed in 
Ref. [1]) at the surface of the superconductor connected 
to an electron reservoir , i.e., when rie^' = n^^^' = np in 
the superconductor. Here Ra ~ 1/Ga, Rt = 1/Gt- 



spond to normal layers surrounding the superconduct- 
ing island. Information that the island is superconduct- 
ing is included in A± definition through G±: 



,(2) 



Gf' 



(15) 
(16) 



A circuit illustration of this boundary condition is given 
in Figl5^. Electron and hole currents entering the left 
side of the pyramid flow in one normal layer, the right 
currents flow in the other normal layer. —A- resis- 
tance describes Andreev reflection, A^ — quasiparticle 
normal transmission through the superconductor and 
A- is Andreev transmission characteristics. It was as- 
sumed deriving Eas. ([T^ -([T ^ that electrochemical po- 
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tential within the superconductor is constant. This is 
correct if the characteristic distance between NS inter- 
faces of the superconducting island is smaher than the 
electron disbalance characteristic length Aq [S]. It is 
implied that nonequillibrium quasiparticles above the 
gap do not have enough time for energy relaxation at 
their fly through the superconductor, so electrochemi- 
cal potential within the superconductor is constant. Or 
if characteristic bias value SV between adjacent super- 
conducting islands is much smaller than the gap and the 
temperature is much below Tc then most part of the cur- 
rent carry Cooper pairs through superconductors rather 
than quasiparticles above the gap. Then correction to 
the total current in the array from quasiparticle energy 
relaxation in superconductors is expected to be small as 
SV/A <C 1 and discussion related to Xq is not relevant. 

Next important step is expressing all boundary con- 
ditions in terms of the distribution functions ne{h) [1]) 
where 



u„_, u„., u„., u„ u„,, u„^. 



n± 



n± — I dx' 



<^N Jo 



1 



D±ix') 



— 1 I = n± — m±I± 



Here n± — Ue ± Uh and n± = fie ± fih, D±df± = 
I±/o'N- The integral here goes from the NS boundary 
{x = 0) into the depth of the normal metal. Bound- 
ary conditions Eas. (jl3l) - (|12p written in terms of ne(;i) 
will have the same form if one replaces G± by G± — 
G±/(1 + G±TO±). Then 



1 -I- (m 



(2) 
± 



(17) 



where indices 1,2 correspond to two normal layers con- 
tacting with the superconducting island. It is useful 
to work with n^^f^-^ because then the proximity effect 
renormalization of the boundary resistances, m±, is au- 
tomatically taken into account. 

The next step toward current calculation is to draw 
an effective network that describes MAR in the array of 
the junctions using Eas. p^ - (|T^ and evaluate partial 
currents in this network using to Kirchhoff 's laws. Im- 
portant task is to find voltages at the superconducting 
islands. But this can be found easily in several cases. 
For example, when the array consists of equal SNS junc- 
tions, or if most of the voltage drops at normal layers. 
Below these situations will be discussed. More compli- 
cated cases I leave for extended paper. If currents in 
the MAR network are found then electric current can 
be evaluated, e.g., as follows: 



(18) 




L'„, t;„, ty„, u„ 

Fig.4. MAR in a SNSNS array with equal SNS junc- 
tions. The graph shows the effective circuit for quasi- 
particle currents in energy space. The role of voltages 
here play quasiparticle distribution functions. For ex- 
ample n* is the quasiparticle distribution function de- 
pending from E + nV/2\ U„ — npiE + nV/2). Boxes, 
triangles and ovals play the role of effective resistances 
that come from Usadel equations and their boundary 
conditions, Eqs.([T])-((6|. The oval is the resistance of the 
normal layer. The box and the triangle correspond to 
the pyramidal bridge in Fig[3^ and the three-terminal 
device in Fig[3}3. The upper raw of f/'s correspond to 
the first superconductor of the array, the lower raw - to 
the last superconductor. The thick line in the center of 
the graph represents the superconducting island. 



where Ie{h) correspond to the normal layer at the edge of 
the array contacting with the superconducting reservoir 
with zero voltage. 

The format of the letter does not allow to describe 
complicated arrays here, so one of the simplest SNS ar- 
rays, the SNSNS junction (FigH]), will be considered 
below as an example. It will be shown how to construct 
an effective network that helps to describe its transport 
properties. Transport in more complicated arrays, like 
in FiglT] can be described in a similar manner as for 
SNSNS; it will be demonstrated in an extended version 
of this paper. 

Effective MAR network for a SNSNS junction like 
in Figl2] is shown in FigHl The bias between the su- 
percondcutors at the edges of the array is V. Then the 
bias of the central superconductor is V/2; it follows from 
symmetry reasons. Currents J„ and Jn correspond le 
for lines beginning from Um and ending at C/„2 with 
77.2 > ni and —1^ vice-versa. The role of voltages play 
quasiparticle distribution functions. For example is 
the quasiparticle distribution function depending from 
E + nV/2; [/„ = npiE + nV/2). Bars above the distri- 
bution functions are supposed but not shown explicitly 
in the figure. Boxes, triangles and ovals play the role of 
effective resistances (also with bars) that come from Us- 
adel equations and their boundary conditions, Eqs. ((!])- 
The oval is the resistance of the normal layer. The 
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box and the triangle correspond to the pyramidal bridge 
in Fig 12^ and the three-terminal device in FiglSlD. The 
upper raw of f/'s correspond to the first superconductor 
of the array, the lower raw - to the last superconductor. 
The thick line in the center of the graph represents the 
superconducting island. 

Lets find recurrence relations for the currents /„ and 
J„. It follows from Fig |3] that 
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Getting rid of the distribution functions in these set of 
equations one can get: 



In 



1 



a: 



^■"■n-l ^n-1 

— In-2bn-2 



+ {a-n + an-2) 



In+2bn — {U, 

Same equation satisfies J„. Here 

11 1 G 



Un-2 



Gjv 
K = 



r^n 



_L /^n ' 



1 



G" 



2G" 



(19) 

(20) 
(21) 



Ea. (|19p coincides with recurrence relation Eq.(52) from 
Ref. [4 derived for a single SNS junction if I replace in 
Eq.(52) Gn by 



G 



G 



N 



G_ 



N 



2 Gat + G- 



(22) 



It means that SNSNS array behaves similarly as a single 
SNS junction but with energy dependent resistance of 
the normal layer. G^ has singularities at energy corre- 
sponding to the gap edges of the superconducting is- 
land in the center of the SNSNS array. This is the 




Fig. 5. Subharmonic structure in differential conductiv- 
ity of SNSNS junction with the ratio r of SN boundary 
resistance to normal layer resistance Rn, r = 0.2. Peaks 
at half-integer 2/S./V do not appear in dl/dV of SNS 
junctions. 



reason of subharmonic singularities in current voltage 
characteristics if 2A/V — n/2, n = 1,2, .. . instead of 
"usual" positions, 2 A/y — n; FiglS] illustrates it. From 
Eq.(I2D) follows that unusual position of subharmonic 
singularities disappear if Gn — ^ Gn if the resistance 
of the normal layer exceeds the resistance of SN inter- 
faces. Then the central superconducting island of the 
SNSNS array effectively "disappear" . It was checked if 
the exchange field in SFSFS junction splits subharmonic 
structure; there was found no exchange field splitting ef- 
fects because configurations of the ferromagnet and the 
superconductor like in FiglT] does not lead to enough 
exchange field deformation of superconductor DoS near 
SF boundary that is necessary for the splitting effect 
observation. This paragraph is conclusion of the paper. 
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